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subarray  operation.  The  uptime  has  nevertheless  been  95.5%,  which  is 
quite  good  in  consideration  of  the  fact  that  the  operation  is  now  unmanned 
except  for  weekday  daytimes.  Part  of  the  explanation  for  this  is  an 
effective  countermeasure  in  terms  of  an  automatic  telephone  warning 
system  triggered  by  DP  stops.  .The  gradual  transition  towards  unmanned 
operation  has  been  accompanied! by  the  lay-off  of  8 of  the  9 computer 
operators,  who  have  had  little 'problems  in  finding  new  jobs.  The  per- 
formance of  the  data  communications  lines  inside  the  array  as  well  as 
within  the  ARPA  subnetwork  has  been  relatively  satisfactory.  Preliminary 
results  on  the  performance ,6f  the  new  Online  Event  Processor  (which  is 
still  being  developed)  ipdicate  quite  satisfactory  results.  A considerable! 
amount  of  work  has  beer!  invested  in  modifying  and  improving  the  DP.  The 
work  has  partly  been  general  system  changes,  and  in  particular  those 
caused  by  the  reduction  in  array  size  (including  the  development  of  the 
Online  EP) , and  partly  changes  connected  to  the  adaption  of  a new  VELANET  j 
protocol  for  tran-Atlantic  message  exchange.  The  field  work  connected  to  I 
the  closing  down  of  the  15  subarrays  is  now  in  progress,  in  accordance 
with  a detailed  plan.  The  work  continues  with  getting  attenuated  LP  data 
in  parallel  for  all  remaining  subarrays.  Some  modifications  have  been 
made  to  the  analog  SP  station,  and  one  attenuated  SP  channel  has  been 
implemented^^ The  regular  array  maintenance  activity  has  Been  kept  at  its 
usual  level,  and  the  monitoring  schedule  has  been  adapted  to  the  new 
array  configuration.^^. 

The  research  activities'N^re  documented  here  in  eleven  reports  covering  a 
wide  range  of  subjects.  The  first  three  reports  cover  further  investi- 
gations of  the  crustal  and  lithospheric  structure  under  NORSAR,  by  using 
the  spectral  ratio  method  in  studying  Moho  depths  under  NORSAR  (1),  by 
using  relative  amplitude  and  travel  time  anomalies  simultaneously  in  a 
new  and  encouraging  approach  on  investigating  lithospheric  heterogeneities 
(2),  and  by  using  S to  P refracted  waves  in  finding  the  boundary  between 
the  lithosphere  and  the  asthenosphere  (3) . The  following  three  reports 
cover  some  aspects  of  seismicity!  one  report  presents  focal  mechanism 
solutions  for  two  recent  earthquakes  from  Iceland  and  Svalbard  (4) , in 
another  report  the  seismicity  in  the  Svalbard  area  is  further  investi- 
gated, including  an  analysis  of  the  crustal  structure  there  (5) , and  the 
third  seismicity  investigation  is  from  the  Hindu  Kush-Pamir  region, 
where  also  the  lithosphere  has  been  investigated  (6) . Then  follows  a 
report  on  the  P-wave  detectability  of  a large  number  of  WWSSN  stations 
(7) , a new  approach  on  investigating  microseismic  noise  using  auto- 
regressive modelling  (8),  and  an  important  contribution  to  seismic 
discrimination  where  a pattern  recognition  approach  has  been  developed 
(9) . The  following  report  describes  a new  method  for  estimation  of 
spectral  properties  of  spatial  data  (10) , and  the  last  one  is  a short 
note  on  the  problem  of  estimating  the  autocorrelation  function  of  spatio- 
temporal  variables  (11). 
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I.  SUMMARY 

This  report  describes  the  operation  and  research  activities 
at  the  Norwegian  Seismic  Array  (NORSAR)  for  the  period  from 
1 October  1976  to  31  March  1977. 

The  operation  of  the  array  in  this  period  has  been  dominated 
by  the  effects  of  the  change  from  22  to  7 subarrays.  There 
have  been  232  breaks  in  the  operation  of  the  Detection  Processor 
(DP),  121  of  which  have  been  due  to  SPS  malfunctioning.  Moreover, 
many  of  the  breaks  are  tied  to  tests  and  software  changes 
connected  to  the  transition  to  a 7 subarray  operation.  The 
uptime  has  nevertheless  been  95.5%,  which  is  quite  good  in 
consideration  of  the  fact  that  the  operation  is  now  unmanned 
except  for  weekday  daytimes.  Part  of  the  explanation  for  this 
is  an  effective  countermeasure  in  terms  of  an  automatic  tele- 
phone warning  system  triggered  by  DP  stops.  The  gradual  transi- 
tion towards  unmanned  operation  has  been  accompanied  by  the 
lay-off  of  8 of  the  9 computer  operators,  who  have  had  little 
problems  in  finding  new  jobs.  The  performance  of  the  data 
communications  lines  inside  the  array  as  well  as  within  the 
ARPA  subnetwork  has  been  relatively  satisfactory.  Preliminary 
results  on  the  performance  of  the  new  Online  Event  Processor 
(which  is  still  being  developed)  indicate  quite  satisfactory 
results.  A considerable  amount  of  work  has  been  invested 
in  modifying  and  improving  the  DP.  The  work  has  partly  been 
general  system  changes,  and  in  particular  those  caused  by  the 
reduction  in  array  size  (including  the  development  of  the 
Online  EP) , and  partly  changes  connected  to  the  adaption  of 
a new  VELANET  protocol  for  trans-Atlantic  message  exchange. 

The  field  work  connected  to  the  closing  down  of  the  15  subarrays 
is  now  in  progress,  in  accordance  with  a detailed  plan.  The 
work  continues  with  getting  attenuated  LP  data  in  parallel 
for  all  remaining  subarrays.  Some  modifications  have  been 

made  to  the  analog  SP  station,  and  one  attenuated  SP  channel 

* 

has  been  implemented.  The  regular  array  maintenance  activity 
has  been  kept  at  its  usual  level,  and  the  monitoring  schedule 
has  been  adapted  to  the  new  array  configuration. 
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The  research  activities  are  documented  here  in  eleven  reports, 
covering  a wide  range  of  subjects.  The  first  three  reports 
cover  further  investigations  of  the  crustal  and  lithospheric 
structure  under  NORSAR,  by  using  the  spectral  ratio  method  in 
studying  Moho  depths  under  NORSAR  (1),  by  using  relative 
amplitude  and  travel  time  anomalies  simultaneously  in  a new 
and  encouraging  approach  on  investigating  lithospheric 
heterogeneities  (2),  and  by  using  S to  P refracted  waves  in 
finding  the  boundary  between  the  lithosphere  and  the  astheno- 
sphere  (3).  The  following  three  reports  cover  some  aspects 
of  seismicity;  one  report  presents  focal  mechanism  solutions 
for  two  recent  earthquakes  from  Iceland  and  Svalbard  (4) , 
in  another  report  the  seismicity  in  the  Svalbard  area  is 
further  investigated,  including  an  analysis  of  the  crustal 
structure  there  (5),  and  the  third  seismicity  investigation 
is  from  the  Hindu  Kush-Pamir  region,  where  also  the  lithosohere 
has  been  investigated  (6).  Then  follows  a report  on  the  P-wave 
detectability  of  a large  number  of  WWSSN  stations  (7),  a 
new  approach  on  investigating  microseismic  noise  using  auto- 
regressive modelling  (8),  and  an  important  contribution  to 
seismic  discrimination  where  a pattern  recognition  approach 
has  been  developed  (9) . The  following  report  describes  a new 
method  for  estimation  of  spectral  properties  of  spatial  data 
(10) , and  the  last  one  is  a short  note  on  the  problem  of 
estimating  the  autocorrelation  function  of  spatio-temporal 
variables  (11) . 

H.  Bungum 
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II 

II. 1 


OPERATION  OF  ALL  SYSTEMS 

Detection  Processor  (DP)  Operation 
* 

Within  the  reporting  period  there  have  been  232  breaks  in 
the  otherwise  continuous  operation  of  the  Detection  Processor 
(DP)  system.  The  uptime  percentage  is  95.5%,  as  compared 
to  96.6%  for  the  last  reporting  period  ( July-September  1976)  . 

Fig.  II. 1 and  the  accompanying  Table  I I. 1.1  both  show  the 
daily  DP  downtime,  in  hours,  for  the  days  between  1 October 
1976  and  31  March  1977.  The  monthly  recording  times  and  up 
percentages  are  given  in  Table  II. 1.2. 

As  can  be  seen  from  Fig.  II. 1,  there  are  days  with  long  down- 
time intervals  towards  the  end  of  the  period.  In  most  cases, 
these  are  stops  that  have  occurred  during  unmanned  shifts, 
when  the  automatic  telephone  alarm  either  was  not  yet  installed 
or  just  failed.  As  experience  is  gained  with  the  alarm  system, 
we  hope  to  improve  its  performance. 

There  were  two  breaks  on  26  December,  which  were  caused  by  power 
failure  and  Special  Processing  System  (SPS)  malfunctioning. 

These  intervals  together  gave  a total  downtime  of  nearly  14  hours. 

The  232  breaks  occurring  within  the  reporting  period  may  be 
grouped  in  the  following  categories: 


a) 

SPS  malfunctioning 

121 

b) 

Error  on  the  multiplexor  channel 

31 

c) 

Tests 

17 

d) 

Software  related  stops 

15 

e) 

Operation  stops 

13 

f) 

Maintenance 

7 

g) 

CPU  disabled 

7 

h) 

Tape  drive  problems 

6 

i) 

Power  j limps  and  breaks 

5 

j) 

Other  hardware  (TOD,  Air  condition) 

4 

k) 

EOC  unit  problems 

3 

1) 

Unknown 

3 

4 


I 


F ' 


i 


It  is  evident  that  the  source  of  the  major  part  of  the  breaks 
in  the  recording  is  the  SPS . Since  diagnostic  aids  are  very 
limited  for  this  system,  it  is  difficult  to  track  down  any 
intermittent  failures  such  as  those  occurring.  (For  futher 
information  about  the  status  of  the  SPS,  see  elsev/here  in 
this  report.) 

The  error  on  the  multiplexor  channel,  which  occurred  in 
abundance  in  the  beginning  of  this  reporting  period,  has  not 
occurred  after  the  beginning  of  January.  This  is  because  the 
load  on  the  multiplexor  channel  has  been  decreased,  by  moving 
the  output  for  the  printer  to  a tape  on  a selector  channel 
instead,  and  also  by  reserving  the  multiplexor  channel  ex- 
clusively for  ARPANET  message  traffic  during  unmanned  shifts 
(as  described  under  Improvements  and  Modifications) . Thus  we 
have  virtually  eliminated  the  occurrence  of  this  error. 


The  rather  high  number  of  tests  and  software  related  stops 
(program  changes  and  program  errors)  must  be  seen  in  relation 
to  the  system  being  in  a state  of  transition  throughout  this 
period,  with  frequent  modifications  and  improvements.  Also, 

I 

the  technique  of  testing  out  the  new  versions  of  the  system 
* as  the  Secondary  Online  system  (on  the  B-computer)  before 

putting  it  into  operation  has  been  used  very  little,  because 
of  the  risk  of  bringing  down  the  SPS  while  initiating  the 
Secondary  Online  system. 


The  category  g)  (CPU  disabled)  is  almost  certainly  a hardware 
error,  occurring  while  executing  one  of  the  microcodes  on  the 
A-computer.  At  least  one  of  these  cases  was  positively  identi- 
fied as  such.  All  these  stops  have  occurred  within  the  last 
month,  and  further  investigations  will  probably  give  more 
conclusive  evidence. 

The  total  downtime  for  this  period  was  204  hours  3 minutes. 

The  mean-time-between-failures  was  0.8  days  as  comoared  with 
0.4  days  for  July-September  1976. 


Detection  Processor  Downtime  Oct  1976  - Mar  1977. 
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LIST  OF  bREAKS  IN  DP  PROCESSING  THE:  LAST  HALF-YtAR 
LAY  START  STOP  COMMENTS.... 


302 

1 

6 

1 

21 

SPS 

3C2 

10 

2 

10 

13 

1052  MAINTENANCE 

3C3 

0 

52 

1 

6 

SPS 

3Gb 

9 

U, 

■A 

9 

2 9 

MOVE  PRINTOUT  TO  1403 

308 

14 

1 1 

14 

16 

MOVE  PRINTOUT  TO  TAPE 

312 

17 

0 

17 

15 

SPS 

313 

0 

0 

34 

SPS 

313 

0 

3 9 

10 

10 

SPS  AND  CHANGE  Trie  1052 

313 

2 3 

3 7 

24 

0 

SPS 

314 

0 

0 

0 

b 

SPS 

314 

3 

4 

3 

13 

SPS 

3 14 

13 

12 

’ 3 

30 

SPS 

315 

5 

6 

3 

2 7 

PROG  CHECK. 

3x& 

2 

43 

3 

7 

SPS 

316 

18 

A 

lb 

21 

SPS 

32  4 

11 

K 

1 1 

1 5 

MOVE  PRINTOUT  TO  1403 

324 

11 

5 3 

12 

8 

MOVE  PRINTOUT  TU  TAPE 

329 

12 

•+  ■> 

13 

2 

PROGRAM  Tc ST 

330 

7 

4 3 

8 

2 

PROGRAM  TEST 

330 

13 

i.  0 

• j 

A - 

a.  \j 

PROGRAM  TEST 

331 

13 

Q 

13 

10 

PROGRAM  TEST 

3 3 4 

10 

3 b 

10 

4 7 

OPERATOR  ERROR 

335 

8 

0 

8 

13 

SPS 

336 

10 

lb 

10 

33 

POWER  JUMP 

3 36 

14 

3 9 

14 

45 

LUC  POWER  ON 

3 3 6 

22 

44 

24 

0 

MPX7LATL  (UNMANNED) 

33  7 

t \ 
0 

0 

6 

58 

MPX/LATt  (UNMANNED) 

337 

14 

14 

19 

PROGRAM  CHECK 

337 

23 

34 

24 

0 

MPX/LATt  (UNMANNED) 

33b 

0 

0 

0 

52 

MPX/LATE  (UNMANNED) 

343 

8 

52 

9 

3 

PLOTTER  MAINTENANCE 

3^3 

1 3 

c 7 

15 

39 

SPS 

34  5 

8 

24 

e 

30 

EOC  DOWN 

}4« 

14 

37 

14 

42 

SPS  maintenance 

350 

10 

12 

10 

30 

PROG  CHECK 

3 30 

1 <1 

2 7 

1 c. 

4d 

PRUGRAM  CHANGE 

3 50 

16 

48 

17 

5 

PROGRAM  CHANGE 

j 33 

6 

16 

15 

IS 

POWER  SUPPLY  MAI  NT. 

353 

16 

21 

16 

29 

cOC  POWER  ON 

334 

T 

A. 

3 

7 

14 

MPX/LATE 

334 

15 

4 6 

16 

0 

MPx/LATE 

2 34 

23 

4 f; 

<.4 

n 

MPX/LATE 

355 

0 

0 

6 

22 

MPX/LATE 

3 pb 

11 

36 

12 

L 6 

PROGRAM  CHANGE 

337 

1 3 

30 

13 

38 

TOD  CORRECTION 

361 

1 

4 0 

10 

24 

POWER  QP.lAK 

361 

13 

24 

18 

24 

SPS 

56  6 

18 

i 6 

j.  4 

0 

SPS 
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LIST  OF  bREAKS  ['1  DP  PROCESSING  1 HE  LAST  HALF-YEAR 


LAY 

START 

S TOP 

COMMENTS. 

1 

0 

G 

0 

1 3 

SPS 

4. 

G 

t X 

x 

2 

' j r*  'i 

*. 

* L 

X.  > 

11 

23 

TIMING  INCORRECT 

) 

X 

i O 

J 

LI 

2 3 

timing  incorrect 

2 

8 

.'i  r 

?> 

43 

IIMiNG  INCORRECT 

3 

9 

L -> 

9 

24 

PRU-jKAM  CnANGt 

1C 

A 2 

1 A 

X u 

4 7 

jP  s 

u 

Ih 

•*7  I 

’ *♦. 

61 

MP  K/LA  Tt 

4 

IS 

X / 

X L,’ 

. 4 

SPS 

•V 

80 

i.  «. 

X G 

6 0 

3P  s 

t, 

0 

v-  f' 

6 

8 

IPS  ( UNMANNED > 

VP 

21 

2 / 

21 

3 7 

SPS  (ALARM) 

9 

16 

<4  _ 

1 6 

4 8 

SPS 

9 

22 

4.  A 

82 

4 6 

>PS  (ALARM) 

*0 

! 

A) 

( 

S4- 

ii1  j 

1 0 

r) 

a 

i G 

S3 

St' s 

iu 

C V. J 

> L> 

<•0 

39 

PROG  CHECK 

x x 

x 

f 

SO 

‘•*1 1 W 8 ft  Bk  x.  AK 

X X 

b 

0 

3 6 

SPS 

‘ T 
L X 

13 

> A 

x J 

i>l 

PF  J jR  AM  CHA  ILL 

*.  A 

17 

IS 

14 

TAPE  DRIVE 

1 <- 

x X 

4.  /T 

X x 

2 7 

UNKNOWN  REASON 

i 3 

10 

6 

10 

20 

SPS 

x "J 

8 

*T  ‘ . 

b 

SPS  (ALARM) 

i.  7 

0 

A t 

4. 

Q 

SPS  (ALARM) 

19 

t * 

: 9 

5C 

TAi't  DPIVl  CONTROLLER 

i G 

r*t 

£ u 

t 

43 

SPS 

*0 

* ; 

a.  w 

lb 

»PS  MA  l '1  i i_  NANCE- 

X G 

1 2 

T 

12 

’ "> 

«.  X 

SP  S 

cO 

12 

4 

».  4 

. 7 

3P  3 

c c 

21 

*' 

21 

3 7 

SPS 

2 / 

22 

5 . 

• . *L 

S3 

SPS 

c 3 

n 

- 3 

30 

SPS 

c.  i 

6 

/ 

6 

44 

SPS 

6.  3 

» 

: o 

8 

23 

SPS 

/ 1 

8 

«-  y 

9 

16 

SPS 

c 3 

10 

X G 

10 

37 

SPS 

X J 

1 0 

3 g 

I G 

34 

s»s 

8<* 

14 

X < * 

G 

3,; 

jH  S 

c4 

16 

19 

1 G 

S J 

SPS 

84 

23 

6 

2 3 

ia 

si's 

2 G 

n 

*.  ^ 

9 

4A 

SI'S 
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( LIST  OF  bRfcAKS  IN  DP  PROCESSING  THE  LAST  HALF- YfcAR 

faAY  START  STUP  COMMENTS* 


f 


A 


Zb 

13 

13 

13 

33 

SPS 

Zb 

13 

A fa 

14 

L 

SPS 

27 

fa 

1 1 

7 

38 

AIR  CONDITION  FAILURE 

2 7 

b 

3 3 

9 

5 

TEMPERATURE  HIGH 

28 

17 

A 8 

18 

fa 

UNKNOWN  REASON 

26 

lb 

43 

18 

51 

UNKNOWN  REASON 

29 

13 

a: 

1A 

IE 

SPS 

30 

A 

1 

A 

2 ^ 

SPS 

30 

7 

24 

7 

51 

SPS 

33 

17 

4 3 

17 

A9 

SPS  ROS  WORD 

3 3 

19 

2 

20 

15 

SPS  ROS  WORD 

3b 

A 

6 

7 

A3 

SPS  (ALARM  FAILURE) 

36 

8 

lb 

6 

36 

SPS 

36 

21 

'3 

2A 

0 

SPS  (ALARM  FAILURL) 

3 7 

0 

0 

8 

3fa 

SPS  (ALARM  FAILURE) 

3 7 

17 

A 2 

17 

AO 

SPS 

38 

12 

58 

13 

3 

SPS 

39 

9 

^ ? 
A *. 

9 

23 

SPS 

AO 

6 

1 

fa 

12 

SPS 

aO 

1 1 

fa 

11 

11 

SPS 

AO 

1A 

A4 

1A 

37 

SPS 

AO 

21 

A L 

22 

A 

SPS 

A3 

1 

A A 

1 

57 

SPS 

A3 

lfa 

lfa 

: fa 

51 

SPS 

45 

2 

31 

2 

50 

SPS 

A3 

9 

i.4 

9 

20 

SPS 

Afa 

9 

A4 

11 

28 

SPS 

A b 

1A 

2 7 

1 t> 

21 

SPS 

45 

lfa 

1 

lfa 

9 

SPS 

A 5 

19 

2 

19 

15 

SPS 

4fa 

20 

39 

21 

3 A 

SPS 

4fa 

23 

2 7 

23 

Afa 

SPS 

A6 

fc 

2 8 

9 

A3 

SPS 

A6 

lfa 

Afa 

1 fa 

54 

SPS 

A6 

1 7 

38 

17 

49 

SPS 

A6 

20 

3u 

20 

36 

SPS 

A7 

17 

3 A 

17 

Ai 

SPS 

A9 

3 

3 3 

A 

13 

SPS 

A9 

7 

X 

7 

10 

SPS 

A 9 

9 

A . 

9 

26 

SPS 

A9 

li 

24- 

11 

28 

SPS 

5 1 

11 

fa 

L x 

12 

SPS 

3a 

fa 

J fa 

t. 

58 

TAPE  DRIVE 

3 3 

fa 

f 

t> 

27 

TAPE  DRIVE 

33 

17 

50 

17 

55 

SPS 
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lIa!  UF  UREAI-S  IN  UP  PROCESSING  THE  LAST  HALF-YtAR 


. M Y 

b T mK 

T 

stop 

COMMENTS. 

:>4 

i j 

i i 

X .2 

23 

SPS 

lb 

j * 

2-  2 

3 ') 

— L. 

SPS  (ALARM  FAIL  UKr.  ) 

22 

5 3 

0 

SPS 

3A 

23 

3 7 

- A 

0 

SPS  (ALARM  FAILURE) 

'J  J 

0 

c 

0 

A 

SR  S (ALARM  FAILURE  ) 

33 

17 

^ 4 

17 

4 J> 

SR  S 

.jb 

12 

' A 

i z 

*»6 

SP  s 

36 

1A 

A 3 

1 4 

5 7 

SPS 

3 8 

5 

4 4 

s 

' U 
_ ^ 

SPo 

39 

A 

: j 

4 

A0 

SPS 

> ) 

1 

./  L 

1 0 

17 

:>PS 

j ) 

10 

2 * 

10 

51 

SPS 

j ) 

- £ 

f> 

X 2_ 

36 

>P  s 

j9 

13 

C.  4. 

^ J 

26 

SP  s 

oC 

vJ 

t 0 

0 

A j 

SPS 

ou 

■7 

J _ 

9 

3 9 

SPS 

c 0 

• 2- 

2 . 

4.  j 

26 

PROGRAM  TlST 

oo 

1A 

*5 

.A 

28 

PROGRAM  TlST 

oO 

1A 

4 6 

„ J 

3 

PRO. .RAM  TEST 

6 J 

15 

3 

«.  -> 

a5 

pro;, pam  (lst 
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1 - 
i.  2 

> 4 

17 

V r 

PROGRAM  ft.  ST 

61 

0 

4 * 

9 

A 9 

PROGRAM  TlST 
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'/. 

1 7 
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1 7 

A 
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6 
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7 
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3 / 

4 'J 
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o3 

1 3 
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- ^ 

c 

PROGRAM  T'-ST 

O 5 

*.  J 

<_  * 

. 7 

59 

PROGRAM  I ST 

63 

1 8 

.< 

_ 7 

22 

PROGRAM  TtST 

L tj 

▲ <. 

_>  3 

1 3 

_ o 

< * rM  J hi  x/L  L \ 

L 7 

1 "j 

5 7 

1 

i » ^ 

i ri 
1 X/ 

i.pij  uiSAbLro 
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i -V 
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i 1 

SP  . 

o ) 

4 

L. 

7 
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) 

->  / 
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2 2 
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2 7 
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•» 

4 d 

>p  J 

tz 

8 

Z ) 

A. 
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1 c 

2 0 

} 
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1 *J 

,PS 

!<. 

20 

4 .■ 

X 1 

2 A 

SPS 

7 3 

2— 

* 

A3 

SP  j 

7 3 

lA 

. 7 

£, 

^ » 

w X 
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7a 

11 

4 „ 

* i. 

16 

SPS 

7 3 

i d 

* > 

- d 
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LIST  U F BREAKS  IN  UP  PROCESSING  THE  LAST  HALF-YEAR 


DAY 

START 

STOP 

COMMENTS. 

77 

16 

20 

17 

1 ^ 

s 

76 

9 

9 

A9 

SPS 

79 

15 

3 5 

15 

56 

SPS 

79 

16 

i 3 

16 

25 

SPS 

79 

17 

Ao 

2 2 

16 

SPS  maintenance 

«0 

3 

AC 

A 

7 

SPS 

d 3 

7 

11 

6 

2 

SPS 

. 63 

15 

23 

1 3 

3A 

SPS 

as 

20 

7 

2 A 

0 

SPS  (ALARM  FAILURE) 

66 

0 

0 

6 

53 

SPS  (ALARM  FAILURE) 

c 7 

11 

.6 

1 1 

■A. 

3 9 

CPU  DISABLED 

o 7 

13 

36 

13 

A 6 

POWER  JUMP 

1.  d9 

fc 

'♦  J 

6 

A 9 

TAPI  DRIVES 

: J 

By  J 
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II . 2 Event  Processor  (EP)  Operation 

The  regular  operation  and  maintenance  of  this  system  was 
discontinued  after  1 October.  However,  some  processing  has 
still  been  done  as  can  be  seen  from  Table  II. 1.2. 


D.  Rieber-Mohn 

II . 3 NORSAR  Data  Processing  Center  (NDPC)  Operation 
Data_Center 

The  reduction  in  the  number  of  computer  operators  from  9 to  1 
took  place  in  the  period,  accompanied  by  a gradual  change 
from  manned  to  unmanned  operation.  Up  to  20  November  a full 
24-hour  shift  plan  was  used  with  lack  of  operators  for  only  a 
few  shifts.  From  20  November  to  15  January  there  was  manned 
operation  in  the  period  0700  to  2200  throughout  the  week. 

The  rest  of  the  period  there  was  manned  operation  from  0700  to 
2200  Monday  through  Friday  and  7 hours  on  Saturday  and 
Sunday.  As  for  the  operators  who  had  to  leave,  five  were 
engaged  as  operators  and  one  as  a programmer  at  other  data 
processing  centers  in  the  area. 

In  order  to  be  able  to  restart  the  Detection  Processor  when  it 
stops  due  to  some  failure  in  unmanned  periods,  an  automatic 
calling  device  has  been  installed.  When  DP  stops,  the  calling 
device  is  initiated  and  a telephone  call  is  made  to  the  person 
whose  telephone  number  has  been  programmed  into  the  calling 
device.  This,  together  with  some  changes  in  the  DP  program, 
is  the  reason  why  DP  uptime  is  as  high  as  95.5%. 

Maintenance  of  computer  equipment  is  in  general  still 
taken  care  of  by  subcontractors.  Project  personnel 
have,  however,  extended  the  area  of  participation  and 
cover,  in  addition  to  the  SPS  and  Exnerinental  Opera- 
tions Console  (EOC) , also  Tape  Units.  Our  aim  is  to  reduce 
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the  subcontractors'  maintenance  engagement,  not  only  to  special 
equipment  like  SPS,  EOC , etc.,  but  also  tape  unis,  certain 
controllers  and  equipment  with  a high  degree  of  reliability. 

Data_Communications 

The  number  of  subarrays  was  reduced  from  22  to  7 on  30 
September  1976  and  the  number  of  communications  circuits  was 
reduced  accordingly.  The  performance  of  the  remaining  circuits 
has  been  relatively  satisfactory,  even  though  groups  of 
subarrays  have  been  affected  by  irregularities  in  carrier 
frequency  equipment,  power  failures,  etc. 

Single  subarray  communication  systems  have  been  affected  for 
other  reasons,  such  as  broken  telephone  cables,  power  failures, 
faulty  equalizers,  and  heavy  damping  on  special  frequencies. 
These  incidents  normally  caused  longer  outages.  Subarrays 
which  have  been  particularly  affected  are: 


OlA 

Week 

2/77 

15%  outage 

OlA 

fl 

6/77 

14%  outage 

02B 

II 

2/77 

25.7%  outage 

02C 

II 

42,44/76 

20.8  and  37.4%  (reduced  performance, 

< 200  errors  per  16  2/3  minutes) 

02C 

If 

51/76 

22.8% 

0 2C 

II 

2/77 

5.1%  outage 

02C 

II 

5,6/77 

7.2,  8.6%  (reduced  performance, 

< 200  errors  per  16  2/3  minutes) 

The  performance  of  the  modems  and  other  communications  hardware 
has  been  satisfactory.  In  March  1977  the  Central  Terminal 
Vault  (CTV)  modem  at  06C  failed  due  to  a faulty  card.  In 
order  to  improve  the  reliability  of  the  existing  communications 
system,  the  Norwegian  Telegraph  Administration  (NTA)  has 
remeasured  all  circuits  and  new  equalizers/amplifiers  will 

be  installed. 

Table  II. 3.1  indicates  degraded  performance  and  outages  of 
the  different  subarray  communications  circuits. 
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The_ARPA_Subnetwork_^TIP_ to_TIP , i • e . x TIP-modem/com_.  lines/ 
modem-TIP^ 

The  London  Circuit 

The  operation  has  been  satisfactory  apart  from  a few  incidents 
with  'Carrier  Loss',  indications. 

The  SDAC  Circuit 

In  October  1976  deteriorated  data  were  reported  relatively 
often.  Tests  and  checks  were  run  with  participation  of  such 
agencies  as  NCC,  ITT,  and  NTA.  The  rest  of  the  reporting 
period  was  far  more  reliable,  although  'Carrier  Loss'  and 
flashing  in  'Marginal  Circuit'  indicators  (both  are  modem 
indicators)  were  observed  from  time  to  time. 


The  Terminal  Interface  Message  Processor  (TIP)  (Including 
Modems  and  Interfaces) 

* The  TIP  equipment  operation  has  been  satisfactory,  apart  from  a 

period  in  October  1976  with  particularly  long  outages. 

, These  outages  were  mainly  caused  by  a Modem  Interface  (no. 

1)  where  a defective  card  had  to  be  replaced.  In  addition 
the  CPU  and  the  power  supply  cooling  fans  were  out  of  operation, 
the  former  due  to  a blown  fuse  and  the  latter  due  to  a 
burned  field  winding.  Lack  of  cooling  was  the  probable 
reason  behind  the  Modem  Interface  trouble,  which  occurred  in 
the  period  22-25  October  when  the  TIP,  apart  from  relatively 
• short  periods,  was  down  most  of  the  time.  Also  26  October 

was  troublesome  with  respect  to  the  TIP.  There  were  several 
reasons  for  that:  TIP  crashed  due  to  lost  'key-in-loader' 
(software) . Also  start  button  failed  causing  trouble  when 
trying  to  read  and  restore  ' key-in- loader ' . In  addition, 

Modem  Interface  No.  1 operation  was  intermittent  and  so  was 
the  SDAC  communications  line. 
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In  December  the  TIP  was  turned  off  for  10  hours  due  to 
installation  of  a special  'Relay  Supervising  Unit'  in  the  AC 
mains  to  the  computer  hall.  At  the  same  time  the  TIP  power 
cable  was  replaced  in  order  to  protect  the  TIP  system  against 
longer  power  outages  during  unmanned  periods  (when  'Supervising 
Relay'  operates  and  has  to  be  reset  manually) , and  because  of 
the  need  for  a power  cable  able  to  carry  heavier  loads. 

Apart  from  incidents  mentioned  above,  and  minor  trouble 
4 November  1976  (due  to  software  and  a blown  fuse  in  the 
SDAC  modem),  the  TIP  has  performed  well.  The  TIP  has,  of 
course,  stopped  later,  but  restart  could  be  done  in  a straight- 
forward manner  most  of  the  times.  Average  number  of  outages 
per  month  in  the  reporting  period  is  3.3. 

J.  Torstveit 
O.A.  Hansen 

; II. 4 ARPANET 

The  Norwegian  Defence  Research  Establishment  (NDRE)  has  attached 
t its  system  to  the  Very-Distant-Host  ( VDH)  port  during  this 

period,  and  performed  some  measurements  on  this  connection. 

f 

However,  regular  operation  of  this  system  as  a Host  on  the 
network  has  not  yet  started.  Otherwise,  the  attachment  configura- 
•,  tion  and  the  use  of  the  ARPANET  node  has  remained  the  same. 


D.  Rieber-Mohn 
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III.  ARRAY  PERFORMANCE 

The  regular  operation  of  the  old  Event  Processor  was  dis- 
continued on  1 October  1976,  when  the  array  was  reduced  from 
22  to  7 subarrays.  At  the  same  time,  the  production  of  the 
NORSAR  Weekly  Event  Summary  (or  the  seismic  bulletin)  was 
stopped.  The  EP  system,  which  has  demanded  significant  efforts 
both  in  terms  of  computer  capacity  and  personnel  will  be 
replaced  by  an  Online  Event  Processor  (ONEP)  which  is  still 
being  developed.  A preliminary  version  of  the  ONEP  was  imple- 
mented on  24  December  1976.  In  order  to  have  some  data  avail- 
able for  the  largest  events,  and  in  order  to  test  the  performance 
of  the  ONEP,  the  regular  EP  was  operated  from  19  October 
to  31  December  1976,  and  for  some  shorter  intervals  in 
January  1977.  The  Fast  EP  has  also  been  used  at  several 
instances  for  interesting  events,  in  particular  presumed 
nuclear  explosions. 

We  have  conducted  some  analysis  of  the  performance  of  the 
ONEP.  In  a comparison  with  the  old  EP  (for  part  of  December 
1976)  it  was  found  to  perform  quite  favorably.  The  beam  reduc- 
tion procedure  is  somewhat  different  (for  example,  the  ONEP 
does  not  permit  more  than  one  processing  or  effort  at  each 
event)  , and  the  ONEP  was  not  found  to  have  more  cases  in  which 
one  missed  the  epicenter  completely.  The  ONEP,  of  course, 
is  bound  to  choose  a beam  location  for  each  event  (no  refine- 
ment) , but  even  this  is  not  much  of  a problem  since  the  loca- 
tion capability  in  most  regions  using  only  7 subarrays  is  not 
much  better  than  the  average  distance  between  beams.  Using 
a random  selection  of  about  25  events,  it  was  found  that  the 
average  location  difference  between  ONEP  and  EP  was  less  than 
2°  both  in  latitude  and  longitude. 


In  looking  separately  at  the  ONEP  results  for  the  time  period 
up  to  1 March  1977,  it  was  found  that  in  average  10  events 
per  day  could  easily  be  identified  as  either  earthquakes 
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or  explosions.  Between  1 and  2 of  these  were  located  obviously 
incorrectly  and  woul''-  need  to  be  processed  over  again  if 
they  should  be  used,  say,  for  a seismic  bulletin.  For  the 
purpose  of  testing  and  evaluation,  the  threshold  was  kept 
so  low  that  an  average  of  27  events  were  selected  and  plotted 
each  day.  This  could  be  done  easily  because  the  size  of  the 
plots  has  been  reduced  from  19(27)  x 11  inches  in  the  old 
EP  to  5 x 10  inches  in  the  Online  EP . 


H . Bungum 
P.  Engebretsen 
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IMPROVEMENTS  AND  MODIFICATIONS 

Detection  Processor 

The  work  on  the  Detection  Processor  system  within  this  period 

may  be  separated  into  two  (partially  overlapping)  categories: 

a)  general  error  corrections  and  system  improvements,  and 

b)  building  a version  of  the  system  which  is  adapted  to  the 
new  VELANET  protocol  for  message  exchange  via  the  ARPA 
network,  but  which  otherwise  incorporates  the  old  version. 

Under  category  a) , the  following  changes  may  be  listed: 

Two  new  operator  commands  were  added  to  the  system  in 
November.  These  commands  make  it  possible  to  switch 
the  output  intended  for  the  printer  to  a specially 
assigned  tape  ('noprt'),  or  back  to  the  printer  ('print'). 
The  default  option  will  be  that  printer  output  is  written 
to  tape,  since  this  is  more  suitable  for  unmanned 
operations.  Implementation  of  this  set  of  operator 
commands  required  modifications  in  the  code  for  the 
PNOPER2  and  LNMSGR  processors  in  the  Message  Task,  as 
well  as  in  the  CASSGN-routine  of  the  DP-Monitor.  The 
CASSGN  routine  is  now  able  to  change  assignments  for 
non- tape  devices  as  well,  and  to  return  the  old  assign- 
ment to  the  calling  program.  A new  routine,  SWTCHMSG, 
was  introduced  (as  overlay  NDPC2S25)  in  order  to  interpret 
the  new  set  of  commands  and  to  call  the  CASSGN-routine 
accordingly . 

Changes  were  made  to  the  PNRERR  processor  in  November, 
so  that  the  periodic  Error  Report  now  is  printed  for 
valid  (i.e.,  remaining)  subarrays  only,  and  only  if  any 
error  count  for  a valid  subarray  is  nonzero. 

An  improvement  of  the  Host-Going-Down  logic  was  done 
to  the  MONSER  routine  in  the  Network  task.  The  routine 
now  checks  for  ' Host-Going-Down ' signalled  from  the  Data 
Acquisition  Task  on  every  call  (not  only  after  a re- 
activation) , and  sends  out  the  ' Host-Going-Down ' message 
to  the  IMP  before  the  system  goes  down.  


I 


21 


I 

( 


\ 

v 

i 


In  December,  the  first  version  of  the  Online  Event 
Processor  was  developed  and  implemented  into  the 
system.  The  introduction  of  this  new  processor  as  part 
of  the  Detection  Task,  with  its  parameter  lists  and 
buffers,  required  a re-organization  of  the  core  layout. 

A new  phase  was  introduced,  which  contains  the  Tape, 

Detection  and  Network  tasks.  This  phase  is  loaded  during 
system  initialization,  but  after  the  SPS  blocks  have 
been  read  from  the  Core  Image  tape  and  transferred  to 
the  SPS.  The  loading  is  done  by  PNRINT.  Extra  code  was 
added  to  the  PNHRTP  processor,  to  extract  high  rate  data 
from  the  High  Rate  Data  Block  and  insert  into  the  input 
buffer  for  the  Online  Event  Processor.  This  processor, 
called  PNREPX , was  introduced  into  the  Detection  Task, 
and  minor  modif ications  and  adjustments  were  made  to 
other  components  of  this  task  (SNDEPR,  PNRICX,  PNREDS ) . 

A new  Core  Image  tape  was  prepared,  incorporating 
parameters  and  buffers  for  the  PNREPX  processor. 

In  accordance  with  the  general  philosophy  of  frugal 
tape  use,  the  interval  for  writing  LTA-data  to  the 

| 

Detection  Log  tape  was  changed  to  10  minutes  (SNRDAQ) . 

By  the  end  of  December,  two  new  sets  of  operator  commands 
were  implemented.  The  first  set  enables  the  operator 

1 

to  start  ('monit')  and  stop  ('nomon')  logging  of  outgoing 
ARPANET  message  headers  onto  the  Detection  Log  tape,  , 1 

together  with  the  time  of  write  (to  TIP)  completion. 

(The  Detection  Log  tape  may  then  later  be  used  as  input 

for  an  off-line  program  to  give  information  about  the  1 

outgoing  message  flow.) 

The  second  operator  command  set  implemented  ( 1 nomux ' , 

'muxit')  enables  the  operator  to  re-direct  permanently 
all  output  intended  originally  for  the  1052  console 
typewriter  or  the  2540  card  punch,  to  the  printer  instead 
(or  to  the  tape  receiving  printer  output) , or  reversing 
this  option.  The  'nomux'  command  will  be  typed  in  at  the 

* 

A 
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start  of  intervals  of  unattended  operation,  thus 
forcing  all  system  status  output  to  one  (tape)  unit, 
and  leaving  the  multiplexor  channel  free  for  network 
traffic  only.  This  also  reduces  considerably  the  probability 
of  occurrence  of  a 'Multiplexor/Late'  hardware  error 
during  unattended  periods. 

The  new  operator  commands  required  changes  in  the  PN0PER2 
processor,  as  well  as  additions  in  the  Message  Task 
Local  Common  (LNMSGR) . New  subroutines  (MONMSG  and 
MUXMSG)  for  interpreting  the  commands  were  introduced 
as  additional  overlays  (NDPC2$26  and  NDPC2$27) . 

An  error  in  the  PNRICX  processor  in  the  Detection  task 
leading  to  incorrect  end  times  for  incoherent  detections, 
was  discovered  and  corrected  in  January.  Also,  detections 
having  an  LTA  value  equal  to  zero  (occurring  sporadically) , 
were  discarded  ( PNREDS ) . 

- The  Physical  Unit  Block  (PUB)  Table  (AFMPUB)  was  expanded, 
in  January,  to  its  maximum  of  22  entries.  This  allows  a 
pool  of  11  scratch  tapes  to  be  set  up  initially.  Since 
only  3 tapes  are  in  use  at  any  time  (High  Rate,  Low  Rate 
and  Det  Log) , this  allows  for  an  unattended  period  of 
more  than  24  hours. 

Under  category  b) , (buildup  of  a new  version) , the  following 

was  done: 

The  two-pass  processor  in  the  Detection  Task  (PNRSAD) 
that  generates  outgoing  ARPANET  messages  and  queues  them 
to  the  Network  Task,  was  completely  re-designed.  This 
was  done  in  order  to  adhere  to  the  new  VELANET  message 
format  that  will  be  used  in  the  future,  and  also  to 
transmit  low  gain  long  period  data  from  one  subarray, 
and  optionally  20  Hz  short  period  data  from  a requested 
instrument . 
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In  the  Network  (NCP)  Task,  minor  changes  were  done  to 
the  IMPOUT,  IMPWRITE  and  IMPREAD  subtasks  because  of 
the  new  message  lengths  and  formats.  Also,  Command  Control 
Words  (CCWs)  in  NCP  Common  (NCPCOMM)  were  changed 
because  of  this.  In  the  IMPIN  subtask,  the  logic  for 
verifying  an  incoming  VELANET  message  was  changed  to  be 
in  accordance  with  the  new  format.  Also,  extra  logic 
was  inserted  to  check  for  gaps  in  the  sequence  of  incoming 
messages,  and  to  check  for  re-initialization  of  the  SDAC 
CCP  (the  VELANET  header  count  is  then  zero),  and  to 
notify  the  SNNCPR  sequencer  on  such  occasions.  This  latter 
task  was  modified  to  be  able  to  interpret  the  new  intra- 
task signals  and  to  convert  them  to  messages  to  be 
queued  to  the  High  Priority  Message  Task.  Also,  the  process 
of  re-initialization  of  the  complete  task  (and  thereby 
the  ARPANET-CCP  connection)  was  diversified,  so  that 
the  task  may  be  re-initialized  when  desired  without  any 
wait  time.  This  is  contrary  to  the  case  for  the  'Imp 
Going  Down/Imp  Dead'  situation,  where  2 minutes  pass 
before  re-initialization  is  done. 


Also,  a new  mechanism  was  built  into  the  IMPTIM  subtask 
to  prevent  hangups  of  the  outgoing  ARPANET  message  flow 
after,  say,  re-initialization  of  the  TIP  or  Pluribus 
IMP  (this  situation  occurs  sporadically) . The  length  of 
the  time  interval  that  has  passed  since  the  completion 
of  the  last  write  operation  (to  the  TIP)  is  computed 
every  time  this  subtask  is  activated  (each  minute) . If 
this  time  interval  exceeds  a certain  limit  (say,  20 
seconds),  the  SNNCPR  subtask  is  signalled  to  perform  a 
re-initialization  of  the  complete  task,  thus  restarting 
the  write  as  well  as  the  read  operations. 

Additional  logic  has  been  inserted  into  the  Online  Event 
Processor  (PNREPX)  in  order  to,  on  the  basis  of  the 
beam  number  of  the  event  selected,  to  select  and  load  an 
overlay  module,  containing  a segment  of  a table  of  location 


question  are  extracted  and  inserted  into  the  Short  Period 
Variable  File  (SPVF)  section  of  the  output  buffer  from 


this  processor,  along  with  status  parameters  collected 


from  other  parts  of  the  system. 

The  Message  Task  Common  (LNMSGR)  was  expanded  and  modified 
to  incorporate  message  skeletons  for  the  necessary  new 
system  output  messages.  It  was  decided  that  it  is  better 
to  send  back  a message  explicitly  spelling  out  the  generated 
response  of  an  SDAC  data  request  message,  or  alternatively 
send  back  an  error  message,  than  merely  unwittingly  re- 
flecting the  request  message,  as  has  been  done  up  to 
present.  The  SAACMMSG  processor  was  therefore  changed 
accordingly.  The  REFLECT  segment  was  removed,  while 
the  REFORMAT  segment  was  modified  to  give  a more  detailed 
confirmation  message,  or  an  error  message,  indicating 
the  error  in  the  request  message. 

A mechanism  for  saving  the  last  EPX  number  used  by 
the  Online  Event  Processor  across  system-down  gaps  has 
been  implemented,  thus  maintaining  a one-to-one  corre- 
spondence between  EPX  number  and  processed  events  within, 
say,  a year.  This  mechanism  is  implemented  as  follows: 

The  Online  Event  Processor  (PNREPX)  will,  when  finding 


a new  event,  queue  the  EPX  number  of  this  event  to  the 
Disk  Task.  The  Disk  Task  Sequencer  (SNDISR)  recognizes 


this  type  of  queue  block  as  one  containing  an  EPX 
number  and  writes  it  to  a special  location  on  the 
Shared  Disk  Pack.  Whenever  the  Disk  Task  performs 
initialization  after  system  startup,  it  will  read  the 
EPX  number  from  the  location  on  the  Shared  Pack,  add  IP 
and  initialize  the  EPX  counter  in  the  Online  Event 
Processor  (PNREPX)  with  the  result.  In  this  way,  it 
should  be  possible  to  avoid  ambiguity  with  respect  to 
using  EPX  numbers  as  event  identif icators  within  a period 
of,  say,  a year. 


f 
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New  code  for  computing  the  number  of  subarrays  down  was 
inserted  in  the  PNSADP  and  PNDKIR  processors  in  the 
Disk  Task. 

A table  of  location  parameters  (latitude,  longitude, 
distance,  azimuth,  etc.)  has  been  generated  and  solit 
up  into  phase  overlays,  to  be  loaded  from  PNREPX.  There 
is  one  entry  for  each  beam  number  of  each  partition. 

IV. 2 Event  Processor 

This  sytem  has  been  adapted  to  used  new  format  High  Rate 
tapes  as  input. 

D.  Rieber-Mohn 

1 

IV. 3 Array  Instrumentation  and  Facilities 

As  of  1 October  1976  the  communication  lines  to  the  fifteen 
subarrays  to  be  closed  down  were  switched  off  after  prior 
notice  to  the  local  power  companies  and  NTA . Landowners  of 
nine  of  the  CTV  and  LPV  areas  replied  positivelv  to  a reauest 
of  taking  over  the  use  of  the  vaults,  includino  all  respon- 
sibilities. The  last  six  will  then  be  completely  sealed  and 
secured  with  concrete.  Up  to  now  only  the  CTV/LPV  at  11C 
have  been  sealed.  The  status  of  the  array  reduction  program 
regarding  the  CTV/LPVs  is  given  in  the  Table  IV. 3.1.  All  work 
on  the  WHVs  is  remaining  and  is  planned  completed  during  this 
year. 


Transportation  of  the  equipment  from  the  CTV/LPVs  to  NMC  is 
expected  to  take  three  to  four  weeks,  starting  in  the  begin- 
ning of  May. 


Sub- 

array 

Equipment 

demounted 

Equipment 
taken  to  NMC 

1 

Comments 

| 

03B 

X 

X 

1 

04B 

X * 

X 

Except  LPV  equipment 

05B 

X 

06B 

X 

1 

07B 

X 

01C 

X * 

05C 

X * 

X 

07C 

X * 

X 

08C 

X * 

09C 

X 

10C 

X 

11C 

X 

X 

Secured  with  concrete 

12C 

X 

13C 

X 

14C 

X 

* 

Concrete  sealing  remaining. 


Table  IV. 3.1 

Status  of  the  array  reduction  work 


Attenuated  long  period  channels  were  installed  at  subarray  01B 
24  November  1976.  The  output  from  all  three  LP  seismometers 
are  taken  out  in  parallel  with  the  normal  data  channels  on 
the  DBl  attenuator  card,  attenuated  -40  dB  and  fed  through 
separate  Ithaco  amplifiers  through  three  additional  LP  Line 
Termination  Amplifiers  (LTA)  (in  positions  A8  and  A9  in 
SLEM)  to  multiplexer  channels  30,  31  and  32  for  Vertical, 
North-South  and  East-West  components,  respectively.  The 
sampling  rate  is  0.5  Hz,  sampled  on  Random  Data  Addresses  ( RDA) 
13,  14  and  15  (refer  Fig.  IV. 3.1).  The  test  function  routing 
is  given  in  Table  IV. 3. 2,  and  the  MUX  channels  assignments 
in  Table  IV. 3. 3.  The  separate  Ithacos  are  damping  the  seis- 
mometers with  an  average  percentage  of  12,  and  the  damping 
resistances  have  been  increased  accordingly  to  compensate. 
Investigations  have  been  initiated  to  improve  the  attenuation 
circuits,  that  is,  to  decrease  the  Ithaco  gain  and  alter  the 
attenuation  circuit  components  for  better  damping  effects. 


ncluding  attenuated  LP-channels. 


FREQUENCY  HZ 


Fig.  IV. 3. 4 Magnification  of  ground  motion  relative  to  magnification 
at  1.0  Hz  of  NORSAR  Analog  SP  Station,  March  1977. 
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Frequency 

Hz 

Recoder 

Deflection 

mm 

Magnification 

Magnification  Ratio 
relative  to  Magnification 
at  1 Hz 

0.143 

3.0 

153 

0.003 

0.2 

5.0 

500 

0.010 

0.33 

8.5 

2312 

0.046 

0.5 

13.0 

8117 

0.162 

0.625 

16.  5 

16098 

0.322 

0.8 

19.0 

30372 

0.608 

0.95 

19.5 

43956 

0.878 

1.0 

20.0 

49954 

1.000 

1.05 

20.5 

56451 

1.130 

1.25 

20.0 

78053 

1.562 

1.6 

18.5 

118291 

2.368 

2.0 

15.0 

149861 

3.000 

2.5 

12.0 

187327 

3.750 

3.0 

9.5 

213552 

4.275 

4.0 

7.0 

279741 

5.600 

5.0 

5.0 

312211 

6.250 

6.0 

4.0 

359667 

7.200 

7.0 

3.0 

367160 

7.350 

8.0 

2.0 

319704 

6.400 

8.5 

1.5 

270687 

5.419 

9.0 

1.0 

202313 

4.050 

9.5 

0.75 

169062 

3.384 

10.00 

0.5 

124884 

2.500 

G • I • 103 

y = — 

4 7T 2 *f  2 *m 

where  M = the  seismograph  magnification,  X = Damping  0.73,  f = nat . frequency 
1.05  Hz,  X = Cal.  current  400iA,  a = peak-to-peak  deflection  in  mm  on 
the  seismogram,  y = equivalent  ground  motion  in  microns,  m = seismometer 
mass  825  gr,  G = cal.  coil  motor  constant  0.0326  N/A,  f = cal.  frequency. 


Table  IV. 3. 4 


Magnification  of  NORSAR  Analog  Station,  25  March  1977 

(8  Hz  filter) 


MAINTENANCE  ACTIVITY 


A brief  review  of  the  maintenance  accomplished  at  the  subarrays 
by  the  field  technicians  as  a result  of  the  remote  array 
monitoring  and  routine  inspections  is  given.  The  monitoring 
schedule  has  been  changed  to  fit  the  new  array  configuration, 
but  the  test  frequency  on  the  individual  subarrays  is  little 
changed.  Except  for  the  array  reduction  work,  all  maintenance 
referred  to  is  on  the  new  array  configuration. 

Maintenance  Visits 

Fig.  V.l  shows  the  number  of  visits  to  the  subarrays  in 
the  period.  Excluding  maintenance  visits  on  the  communications 
system,  the  subarrays  have  on  the  average  been  visited  3.7 
times.  The  subarrays  in  the  new  array  configuration  have  been 
visited  3.1  times  on  the  average  (excluding  01B:  1.8  times), 
the  large  number  of  visits  to  01B  are  due  to  cable  breakage 
and  installation  of  attenuated  LP  channels. 


No . of 
Visits 


01B  03B  05B  07B  02C  04C  06C  08C  10C  12C  14C 


Fig.  V.l  Number  of  visits  to  NORSAR  subarrays  1 October  1976  - 31  March  1977 
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Table  V.2  gives  the  number  of  accomplished  adjustments  and 
replacements  of  field  equipment  in  the  array  with  the  exception 
of  those  mentioned  in  Table  V.l. 


Preventive  Maintenance  Projects 


The  preventive  maintenance  work  in  the  new  array  is  described 
in  Table  V.l. 


Disclosed  Malfunctions  on  Instrumentation  and  Electronics 


Table  v.l 


Preventive  Maintenance  Work  in  the  Period 
1 October  1976  - 31  March  1977 


Repl.  Ad] 


Characteristic 


Repl.  fid] 
3 

5 


mometer  Damping 


mometer 


I Ampli- 
fier 
|(RA-5, 
Ithaco) 


DCO 


Filter 
Ch.  gain 


LTA 


LE 
SB  Gen 
RSA/ 
ADC 
EPU 


fable  V.2 


rotal  number  of  required  adjustments  and  replacements  in  the  NORSAR  data 
channels  and  SLEM  electronics. 


Action 


Unit 


No.  of  Actions 


Adjustment  of  SP  DC  offset 


LTA 


Adjustment  of 
channel  gain 


Battery  Maintenance 


Power 
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Malfunction  of  Rectifiers,  Power  Loss,  Cable  Breakages 

There  has  been  no  malfunction  of  rectifiers  in  the  period. 
Power  loss  requiring  action  by  the  field  technicians  or  local 
power  company  are  reported  in  Table  V.3. 


Sub- 

Array 

Fault 

Period  of 
Inoperation 

l 

Comments 

02B 

No  power  at  site 

27-29  Nov  76 

No  line  power 

02B 

II  _ 

12-14  Jan  77 

1 

Table  V. 3 
Power  Loss 


One  cable  breakage  (01B03)  has  been  repaired  requiring  13 
days'  work  of  the  field  technicians.  The  period  of  channel 
inoperation  was  from  11  October  - 21  December  1976. 

Conclusion 

Compared  with  previous  periods,  the  performance  of  the  array 
instrumentation  has  been  stable  and  satisfactory  in  this 
period.  Though  adjustment  of  SP  channels  DC  offset  to  a small 
positive  bias  is  part  of  the  routine  preventive  work,  the 
array  average  DC  offset  is  now  -4  millivolts,  but  no  significant 
change  of  the  drift  rate  towards  negative  offset  is  observed. 

The  NORSAR  Analog  SP  station  was  set  in  operation  again 
9 November  after  installation  of  a new  power  supply  system. 

New  versions  of  the  analyzing  programs  were  tested  out  in  the 
first  half  of  the  period  and  catalogued  on  disk  ;( N — 5 ) CHANEVSP 
30  November  and  CHANEVLP,  SACP  and  MISNO  as  of  12  January  1Q77. 

Alf  Kr.  Nilsen 
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VII 
VII  . 1 


SUMMARY  OF  SPECIAL  TECHNICAL  REPORTS/PAPERS  PREPARED 

NQRSAR  Moho  Depths  from  Snectral  Ratio  Analysis  of  Lena 
Period  P-waves 

Using  the  spectral  ratio  method  discussed  by  Berteussen  (1°76), 
we  have  tried  to  estimate  .Moho  depths  under  NORSAR.  In  the 
processing  of  the  data  there  were  several  problems  which  we 
had  to  consider.  From  the  simulation  done  (Berteussen,  1976) 
we  expected  that  a shortening  of  the  window  length  used  in 
the  calculation  of  the  spectra  should  imply  a sort  of  smoothing 
of  the  spectral  ratios.  For  the  real  data  we  did,  however, 
observe  that  in  a number  of  cases  such  a shortening  implied 
an  unsystematic  movement  or  change  in  the  peaks  that  can  bias 
the  results  seriously  if  care  is  not  taken  in  the  interpretation. 
Avoiding  all  cases  where  such  problems  occurred,  we  ended  up 
with  112  spectral  ratio  estimates.  Using  the  size  of  the 
main  spectral  peak,  we  estimated  the  average  crustal  P-velocity 
for  a one-layer  model  to  be  6.6  km/s.  From  the  shape  of  the 
spectral  ratios  we  calculated  depth  to  Moho  only,  as  the 
variation  in  the  ratio  across  the  array  made  it  obvious  that 
it  would  be  impossible  to  find  a more  detailed  model.  The 
Moho  depths  found  were  then  plotted  on  a diagram  of  the  array 
at  the  point  where  the  respective  signals  were  expected  to 
cross  Moho.  The  data  showed  up  to  be  consistent  in  that  we 
could  not  find  drastically  different  depth  values  plotted  close 
to  each  other.  Thus  a signal  coming  from  a certain  direction 
recorded  at  one  instrument  gave  approximately  the  same  Moho 
depth  as  a signal  from  another  direction  crossing  Moho  at 
the  same  location  and  being  recorded  at  another  instrument. 
Smoothing  the  data  and  making  depth  contours  we  then  get 
the  Moho  interface  displayed  on  Fig.  VII. 1.1.  For  more  details 
see  Berteussen  (1977).  This  interface  agrees  well  with  the 
interface  calculated  by  Berteussen  (1975)  which  minimized  the 
short  period  teleseismic  P-wave  travel  time  residuals  measured 
at  NORSAR,  in  that  it  tends  to  be  deep  on  the  western  part 
of  the  array  (down  to  38  km)  and  shallower  (32  km)  on  the 
eastern  and  especially  southeastern  part  of  it.  It  also  agrees 
principally  with  the  structures  found  by  Aki  et  al  ( 1 Q 7 7 ) for 
the  area  under  NORSAR  using  their  block  model. 


K . -A . Berteussen 
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Fig.  VII. 1.1  NORSAR  Moho  depth  contours  estimated  from  spectral 
ratio  analysis  of  long  period  P-waves. 


Modelling  of  P-wave  Amplitude  Anomalies  in  Terms  of 
Heterogeneities  in  the  Lithosphere 

The  novel  block-inversion  scheme  of  Aki,  Chr istof f ersson 
and  Husebye  (1977)  has  proved  very  successful  in  analysis 
of  relative  travel  time  observations  across  an  array  or 
network  of  seismic  stations.  On  the  other  hand,  most  efforts 
aimed  at  explaining  P-wave  amplitude  anomalies  have  only 
been  moderately  successful.  At  NORSAR  due  consideration  has 
been  given  to  this  kind  of  problem  in  view  of  their  importance 
for  seismic  discrimination  and  at  the  same  time  reflecting 
the  need  for  a better  understanding  of  wave  propagation  effects 
of  lithospheric  heterogeneities.  The  first  step  in  this  direction 
was  a comprehensive  mapping  of  all  subarrav  amplitude  anomalies 
for  all  beam  locations  (Berteussen  and  Husebye,  1974)  Berteussen, 
1975) . The  next  step  was  aimed  at  computing  theoretical  ampli- 
tude distributions  across  the  NORSAR  array,  and  then  comparing 
them  with  the  observational  data.  In  this  respect  we  have 
used  a simplified  model  for  the  3-dimensional  structure  beneath 
NORSAR  the  construction  of  which  was  based  on  the  following 
assumptions:  (1)  The  essential  3-D  structures  responsible 
for  the  observed  large-scale  amplitude  and  time  anomalies 
are  located  in  a single  layer  at  a mean  depth  D below  the 
surface,"  (2)  That  P-waves  entering  this  layer  from  below 
may  be  treated  as  plane  waves  with  uniform  amplitude  in  the 
plane  of  the  wave)  (3)  The  scales  of  the  inhomogeneities  are 
such  that  ray  theory  may  be  applied)  (4)  P to  S conversions 
and  anelastic  effects  may  be  ignored.  We  further  assume  that 
so  far  as  the  travel  time  and  amplitude  anomalies  are  con- 
cerned the  postulated  inhomogeneous  layer  may  be  treated  as 
an  equivalent  thin  lens  located  at  the  mean  depth  D of  the 
inhorr.ogeneities . This  thin  lens  model  was  subsequently  con- 
structed by  projecting  travel  time  residuals  for  all  subarravs 
for  about  200  earthquakes  at  various  depths  ranging  from  50 
to  250  km  and  then  averaging.  This  time  residual  model  was 
then  converted  to  the  equivalent  velocity  anomaly  model  and 
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then  Claerbout's  (1976)  approach  was  used  for  numerical  cal- 
culation of  the  wave  field  on  the  free  surface. 

The  results  obtained  so  far  are  very  encouraging,  i.e.,  the 
differences  between  observed  and  calculated  P-wave  amnlitudes 
are  generally  very  small.  Even  more  startling  results  have 
been  obtained  when  reversing  the  above  procedure,  that  is, 
we  have  been  able  to  reproduce  in  great  detail  the  observed 
travel  time  residual  field  from  the  amplitude  anomaly  observa- 
tions via  a modified  version  of  Poisson's  differential  equation. 

A comprehensive  documentation  of  this  work  will  be  completed 
in  the  near  future. 


R.A.W.  Haddon,  Sydney  University 
E.S.  Husebye 
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Long  period  (2-20  s) , longitudinally  polarized  precursors 
to  direct  S have  been  observed  at  the  NORSAR  array  from  large 
earthquakes  at  epicentral  distances  of  70°-82°.  Normal  phases, 
including  multiple  surface  reflections,  can  be  excluded  on 
the  basis  of  slowness  and  travel  time.  The  slowness  of  these 
precursors  differs  enough  from  that  of  the  background  that  the 
arrival  time  can  be  determined  to  within  a second  from  the 
resulting  discontinuity  in  the  particle  motion.  The  determined 
slowness  (approximately  that  of  direct  S)  and  differential 
travel  time  (~28  s relative  to  S)  suggest  the  interpretation 
of  these  arrivals  as  Sg',  the  S-to-P  refraction  is  at  an 
approximately  horizontal  interface  about  235  km  below  the 
Baltic  Shield  2°-3°  from  NORSAR.  The  relative  polarity  of 
Sp  and  S (in  pahse  on  the  horizontal,  out  of  phase  on  the 
vertical)  indicates  a velocity  decrease  with  increasing 
depth  across  the  interface.  We  therefore  suggest  that  the 
discontinuity  is  the  lithosphere-asthenosphere  boundary. 

J.A.  Snoke,  Carnegie  Institute,  Washington 

I.S.  Sacks, 

E.S.  Husebye 
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VI I . 4 Two  New  Earthquake  Focal  Mechanism  Solutions 

| As  a part  of  our  investigations  of  the  seismicity  of  the 

North  Atlantic  Ocean  and  of  the  Greenland  and  Norwegian 
Seas,  focal  mechanism  solutions  have  been  obtained  for  two 
recent  earthquakes,  one  from  Iceland  and  one  from  Svalbard. 

A complete  report  will  be  published  elsewhere  (Bungum,  1977)  . 

i£§I§D5_5§Eii}2H§]S§_2£_i3_JanuarY_19  76 

The  PDE  solution  for  this  earthquake  is  based  on  160  observa- 
tions and  given  as  66.2°N,  16.6°W,  with  origin  time  13.29.19.5 
and  magnitude  6.0.  The  focal  mechanism  solution  presented  in 
Fig.  VII. 4.1  is  based  on  81  first  motion  readings  from 
predominantly  long  period  and  broad  band  seismograms.  The 
present  focal  mechanism  solution  is  in  general  agreement  with 
i'  those  published  for  an  earthquake  from  28  March  1963  (Sykes, 

1967)  and  one  from  5 May  1969  (Conant,  1972) , all  situated 
at  or  close  to  the  Tj0rnes  Fracture  Zone.  We  suggest  that  the 
, actual  plane  of  faulting  for  the  present  earthquake  is  the  one 

striking  at  131°,  the  movement  is  thereby  right- lateral , which 
is  in  accordance  with  the  concept  of  transform  faulting. 

In  comparing  the  strike  angles  for  these  three  earthquakes, 
we  find  a systematic  increase  as  we  move  from  west  to  east 
. along  the  fracture  zone.  Such  deviations,  which  also  have  been 


, 


observed  along  the  Jan  Mayen  Fracture  Zone  (Bungum  and  Husebye, 
1977) , would  be  consistent  with  the  interpretation  of  the 
Tj0rnes  Fracture  zone  as  an  en-echelon  system  of  north-south 
trending  troughs  and  volcanic  chains,  developed  successively 
as  spreading  axes  across  the  fracture  zone  ( Seemundsson , 1974)  . 

Sva  lbar  d_Ea  r e_o  f _1§_ January_ 1976 

The  PDE  solution  for  this  earthquake  is  based  on  128  observations 
given  as  77.9°N,  18.6°E,  with  origin  time  04.46.24.4  and 
magnitude  5.6.  The  focal  mechanism  solution  presented  in 
Fig.  VII. 4.2  is  based  on  61  first  motion  readings  from 


Fig.  VII. 4.1  Focal  mechanism  solution  for  the  Iceland  earthquake  of 

13  January  1976.  Solid  circles  are  compressions,  open  circles 
dilatations,  and  crosses  indicate  stations  with  questionable 
first  motion. 


Fig.  VII. 4. 2 Focal  mechanism  solution  for  the  Svalbard  earthquake  of 

18  January  1976.  Solid  circles  are  compressions,  open  circles 
dilatations,  and  crosses  indicate  stations  with  questionable 
first  motion. 
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predominantly  long  period  or  broad  band  seismograms.  This 
is  the  first  focal  mechanism  solution  published  for  the 
Svalbard  area,  and  it  is  noteworthy  that  a clear  strike- 
slip  solution  has  been  obtained  for  this  intra-plate  earth- 
quake. Little  information  is  available  bearing  on  the  stress 
pattern  in  this  area,  most  likely  it  is  the  accumulated  con- 
tribution from  several  sources  and  orogenic  cycles.  The 
geologic  trend  in  the  area  (notably  late  Caledonian  and 
Devonian  faulting)  is  NNW-SSE,  and  this  is  just  between  the 
directions  of  the  two  nodal  planes  in  Fig.  VII. 4. 2.  We  return 
briefly  to  this  event  in  Section  VII. 5. 


H.  Bungum 
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VII. 5 Seismicity  and  Crustal  Structure  in  the  Svalbard  Area 

The  seismicity  of  the  Svalbard  area,  characterized  by  a 
definite  continental  structure,  is  interesting  on  account  of 
its  proximity  to  the  active  Knipovich  spreading  ridge  and 
also  in  view  of  the  complicated  tectonic  history  of  the  area 
in  question.  We  have  re-examined  the  seismic  activity  here, 
resulting  in  the  seismicity  map  presented  in  Fig.  VII. 5.1. 

The  map  covers  the  time  period  1957-1975,  the  epicenters  are 
from  ISC,  PDE  and  Sykes  (1965)  (and  in  that  order  of  priority 
if  more  tnan  one  of  these  have  reported  the  event) , and  only 
events  with  at  least  6 reporting  stations  are  accepted.  The 
map  also  contains  seven  fault  plane  solutions  from  the  area. 

The  area  covered  in  Fig.  VII. 5.1  includes  the  northern  part 
of  Mohns  Ridge  (MR) , Knipovich  Ridge  (KR) , Spitsbergen  Fracture 
Zone  (SFZ)  and  the  southern  part  of  Nansen  Ridge.  The  seismic 
activity  along  Mohns  Ridge  east  of  the  Greenwich  meridian  fol- 
lows the  ridge  fairly  closely,  and  the  normal  faulting  of 
Event  5 also  confirms  this  typical  behavior  for  an  active 
spreading  axis  (Sykes,  1967).  The  well-organized  earthquake 
pattern  extends  through  the  area  where  the  axis  'bends'  into 
the  Knipovich  Ridge,  which  is  in  accordance  with  the  observed 
continuity  between  the  two  ridges  (Talwani  and  Eldholm, 
in  press) . The  seismicity  along  Knipovich  Ridge  is  more 
uneven  and  dispersed,  and  the  strike-slip  solution  of  Event 
1 may  suggest  movements  along  a very  local  fracture  zone. 

The  earthquake  activity  around  the  Spitsbergen  Fracture  Zone 
is  also  quite  dispersed,  which  has  been  used  as  indication 
for  an  en-echelon  structure  of  this  transform  fault  (Husebye 
et  al,  1975).  It  is  noticeable  that  the  3 events  for  which 
focal  mechanism  solutions  are  available  (Events  2-4  in  Fig. 

VII. 5.1),  have  strike  directions  coincident  with  those  of 
the  fracture  zone. 


The  intraplate  seismicity  in  Fig.  VII. 5.1  is  confined  to  two 
areas,  the  north-east  coast  of  Greenland  (with  the  normal 
faulting  of  Event  6)  and  the  Svalbard  region.  It  is  seen  from 
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Fig.  VII. 5.1  Earthquake  occurrence  in  the  Greenland/Svalbard  region. 

Only  epicentral  solutions  based  on  at  least  6 stations 

are  used,  and  the  larger  symbols  indicate  a magnitude 

of  at  least  5.0.  Black  and  white  areas  for  the  seven 

focal  mechanism  solutions  indicate  areas  of  compression 

and  dilatation,  respectively.  The  structural  information 

on  Mohns  Ridge  (MR) , Knipovich  Ridge  (KR) , and  the  Spitsbergen 

Fracture  Zone  (SFZ)  is  taken  from  Talwani  and  Eldholm 

(in  press) . 
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Fig.  VII. 5.1  that  the  earthquake  activity  in  the  latter  area 
is  weak,  with  the  exception  of  a cluster  of  events  along  the 
southeastern  coast  of  the  Spitsbergen  island,  where  a NE-SW 
trend  is  observed.  A puzzling  feature  of  the  seismic  activity 
here  is  the  focal  mechanism  solution  for  Event  7,  presented 
in  Section  VII. 4.  The  available  tectonic  information  cannot 
resolve  the  ambiguity  of  the  solution,  as  the  main  trend  of 
the  Spitsbergen  geology  is  NNW-SSE  and  thus  just  between  the 
directions  of  the  two  nodal  planes.  However,  the  seismicity 
information  give  some  support  for  suggesting  that  the  focal 
plane  is  the  one  trending  NE-SW,  since  that  direction  coincides 
with  the  lineation  of  the  earthquake  epicenters  in  the  area. 


In  the  second  part  of  this  study  we  have  used  data  from  the 
WWSSN  station  KBS  (see  Fig.  VII. 5.1)  in  an  investigation  of 
the  lithospheric  structure  in  the  Svalbard  area.  To  this  end 
the  spectral  ratio  method  has  been  employed  for  estimating 
local  crustal  thickness  (Phinney,  1964,'  Berteussen,  in 
press) . We  scanned  all  the  KBS  data  for  the  years  1972-73 
and  found  8 events  suitable  for  analysis.  The  observed 
spectral  ratios  for  these  events  (calculated  for  frequencies 
up  to  0.2  Hz)  were  found  to  fit  the  theoretically  expected 
values  best  for  Moho  depths  around  27  km.  In  fact,  there 
were  two  events  with  26,  one  with  27  and  five  with  28  km. 
However,  since  the  spectral  ratio  is  also  dependent  on 
factors  which  are  kept  constant  in  our  calculations  (in 
particular  the  signal  window  length  and  the  velocity-depth 
function) , we  cannot  claim  a precision  any  better  than  +2 
km  in  our  determination  of  the  depth  to  Moho  under  the 
station  KBS  at  Svalbard. 

H.  Bungum 

E.S.  Husebye 

K.-A.  Berteussen 
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It  has  been  recognized  for  several  years  that  the  Hindu  Kush 
region  contains  a zone  of  anomalous  seismicity  which  is  not 


readily  explained  by  any  simple  geometry  of  plate  convergence. 
The  shallow  seismicity  is  most  probably  associated  with  the 
large-scale  thrusting  exhibited  in  the  Himalayan  orogeny,  but 
beneath  this  lies  an  intermediate  zone  of  very  localized 
seismicity,  perhaps  25  km  wide  and  300  km  long,  trending  4C° 
from  the  shallow  zone.  We  present  a delineation  of  this  zone 
by  seismicity  data,  examine  the  regional  P-wave  velocity 
structure  using  2 separate  methods.  The  first  of  these  uses 
the  now  familiar  3-dimensional  velocity  modelling  orocedure 
of  Aki,  Chr istof f ersson  and  Husebye  on  arrival  time  data  from 
a network  of  stations  in  the  Pamir  to  produce  a localized 
picture  of  the  velocity  structure  in  the  anomalous  seismic 
zone.  The  change  in  direction  of  the  seismicity  at  around 
150  km  depth  is  echoed  by  an  abrupt  change  in  the  velocity 
anomaly  lineations  in  the  same  depth  range.  The  projections 
to  depths  of  100  to  200  km  of  arrival  time  residuals  at  a 
group  of  more  widely  spaced  WWSSN  stations  are  used  to  augment 
the  small-scale  3-D  structure  obtained  by  the  first  method 
with  a less  detailed,  but  large  scale,  2-D  model. 


E.S.  Husebye 

P.C.  England,  Cambridge  University 
A. A.  Lukk,  Institute  of  Physics  of  the 
Earth,  Moscow 
L.P.  Vinnik,  -"- 
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P-wave  Detectability  Study  of  Existing  World-Wide 
Seismograph  Stations 

A detailed  knowledge  of  the  event  detection  capabilities  of 
existing  seismograph  stations  is  a prerequisite  in  order  to 
accurately  assess  the  existing  capabilities  for  monitoring 
underground  nuclear  explosions.  Such  information  can  most 
reliably  be  obtained  by  actually  measuring  the  performance 
of  each  station  over  an  extended  period  of  time,  in  terms 
of  the  number  of  reported  teleseismic  events.  The  method 
of  Ringdal  (1975)  is  well  suited  to  this  kind  of  analysis, 
provided  an  adequate  'reference  system'  (i.e.,  an  organization 
reporting  a comprehensive  event  catalogue  and  associated  magni- 
tudes) is  available.  In  this  respect,  an  analysis  of  the  event 
detection  capabilities  of  seismograph  stations  in  Fennoscandia 
was  undertaken  by  Pirhonen  et  al  (1976),  based  upon  ISC  reported 
events  during  the  years  1964-69.  Fig.  VII. 7.1  shows  the  type 
of  results  obtained  by  this  study,  and  the  excellent  fit  of 
the  detection  curves  shows  that  the  resulting  detectability 
estimates  must  be  considered  very  reliable. 

The  objective  of  the  present  study,  which  is  as  yet  in  its 
initial  stage,  is  to  perform  a comprehensive  evaluation  of 
world-wide  P-wave  detectability  of  existing  seismograph 
stations  in  terms  of  ISC  magnitude.  The  incremental  probability 
G of  detection,  given  event  magnitude  m,  is  assumed  to  be 
of  the  form  (Ringdal,  1975) 

m . 

G(m;u,o)  = / — • exp (- (x-p ) 2/ ( 2o 2 ) ) dx  (1) 

/2tt  • a 

where  the  two  parameters  p and  a characterize  the  detection 
capability.  Assuming,  for  a given  station,  that  the  probability 
of  'false  alarm'  is  equal  to  V„,  and  the  probability  of  the 


DETECTED 


DETECTED 


NOT  DETECTED 


NOT  DETECTED 


<4.0  5.0 

ISC  MAGNITUDE 


MAX  LIKELIHOOD 
CURVE 

90  PER  CENT 
CONFIDENCE  LIMITS 


MAX  LIKELIHOOD 
CURVE 

90  PER  CENT 
CONFIDENCE  LIMITS 
OBSERVED 
DETECTION 
PERCENTAGES 


OBSERVED 

DETECTION 

PERCENTAGES 


ISC  MAGNITUDE 


ISC  MAGNITUDE 


DETECTED 


DETECTED 


■ NOT  DETECTED 


NOT  DETECTED 


t.l.v.l  u, 


5.0 

isc  mnnuDE 


ISC  MAGNITUDE 


MAX  LIKELIHOOD 
CURVE 

90  PER  CENT 
CONFIDENCE  LIMITS 


MAX  cl KEL I HOOD 
CURVE 

90  PER  CENT 
CONFIDENCE  LIMITS 


» 08SERVED 
DETECTION 
PERCENTAGES 


OBSERVED 

DETECTION 

PERCENTAGES 


4.0  5.0  6.0 

ISC  WiNITUDE 


ISC  MAGN I TUD£ 


Teleseismic  (30  -90  ) P-wave  detection  statistics  for 
selected  stations  based  on  the  years  1964-69.  The  upper 
half  is  a histogram  showing  the  reference  event  set  and 
the  number  of  events  actually  detected  for  each  magnitude 
the  lower  half  shows  the  maximum  likelihood  detectability 
curve  and  its  confidence  limits.  The  actual  percentage  of 
detected  events  at  each  magnitude  is  also  shown. 

A.  For  station  NUR  (Nurmi jarvi,  Finland 

B.  For  station  UPP  (Uppsala,  Sweden) 

C.  For  LHN  (Lillehammer , Norway) 

D.  For  COP  (Copenhagen,  Denmark) 
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station  being  out  of  operation  is  equal  to  Vp , it  follows  that 
the  probability  of  finding  a station  detection  that  can  be 
associated  to  a given  event  of  ISC  magnitude  m is: 

P(m)=P(m;u,a,VF,VD)=(l-VD) • [Vp+(1-VF) *G(m;p,o) ] (2) 

For  a given  set  S of  reference  events,  denote  by  D the  set 
of  indices  for  which  the  station  actually  has  a detection. 

The  logarithm  of  the  associated  likelihood  function  then 
becomes  (Ringdal,  1975) 

log  L ( y , a , V , V ) = l P (m . ) + Y (1-P(m,))  (3) 

i£D  1 if D 1 

This  function  may  subsequently  be  maximized  to  obtain  maximum 
likelihood  estimates  of  the  parameters  y , a , V„ , Vr, . In  practice, 

Vp  may  reasonably  be  set  to  0,  since  the  chance  of  a 'random 
association'  in  the  ISC  bulletin  can  be  considered  very  small. 
However,  the  'downtime'  parameter  Vp  cannot  in  general  be 
ignored,  and  must  be  estimated  on  the  basis  of  available  data. 

We  note  in  passing  that  Pirhonen  et  al  (1976)  estimated  VD 
directly  from  operational  statistics  obtained  from  the  ^enno- 
scandian  stations,"  however,  this  procedure  cannot  be  used 
in  the  present  case  as  no  reliable  operational  reports  are 
generally  available. 

Our  initial  detectability  results,  as  presented  here,  have 
been  obtained  by,  for  each  station  reporting  to  the  ISC, 
extracting  those  ISC  reported  events  that  have  been  located 
within  30-90  degrees  epicentral  distance  of  the  station. 

Only  events  with  an  associated  ISC  computed  m^  value  have 
been  used  in  this  reference  set.  Table  1 summarizes  the  observed 
detection  percentages  and  the  estimated  50  and  90  per  cent 
incremental  teleseismic  detection  thresholds  for  a subset 
of  96  stations.  This  subset  consists  of  those  stations  with 
an  estimated  50  per  cent  threshold  of  mj_j=5.0  and  better. 
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For  four  stations  (NORSAR,  LASA,  UBO,  BMO)  no  estimates 
have  been  given  due  to  the  fact  that  these  stations  detect 
essentially  all  ISC  events  in  the  appropriate  distance  range, 

• hence  ISC  is  not  suitable  as  a reference  network  in  these 

cases.  Otherwise  we  note  that  stations  of  high  detectability 
are  found  worldwide,'  e.g.,  Africa  (BNGJ  m,.Q=4.1,  KIC  -4.4), 
South  America  (PNS  - 4.4J  BDF  - 4.6),  Antarctica  (SPA  - 4.0), 
'A§ia  (KBL  - 4.4,'  ELT  - 4.7),  Australia  (ASP  - 4.4,'  WRA  - 4.4), 
besides  several  sensitive  stations  in  North  America  and  Europe. 

Si 

Future  work  on  the  detectability  study  will  concentrate  on 
evaluating  thfe^potential  bias  in  ISC  m^  estimates  and  its 
effects  on  our  detectability  computations.  Also,  a possible 
^ regionalization  of  detection  thresholds  will  be  studied.  More- 

over, an  analysis  will  be\  undertaken  in  order  to  determine 

\ 

to  what  extent  the  detection  of  a given  event  by  stations 
within  a close  distance  of ‘one  another  are  correlated.  This 
is  important  since  global  detectability  studies  (e.g.,  Snell, 
1976)  generally  assume  statistical  independence  in  this 
context . 

I 

F.  Ringdal 

* E.S.  Husebye 

J.  Fyen 

* 
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VII. 8 Autoregressive  Modelling  of  Seismic  Storms 

A nonstationary  autoregressive  (AR)  time  series  model  given 
by 


X(t)  - ax(t)  X(t-l)  - ...  - a (t)  X(t-p)  = 2 ( t ) (1) 

has  been  fitted  to  the  long  period  (LP)  noise  registered  at 
NORSAR.  Here  X(t)  is  the  observed  LP-noise  while  Z(t)  is  a 
white  noise  residual  process  satisfying  E{Z(t)  zTiT }=o^ ( t ) 6 . 

The  LP-noise  is  assumed  to  be  approximately  stationary  within 
time  intervals  of  duration  5 min.  (300  samples),  and  time 
variation  of  the  coefficients  is  introduced  by  fitting  sta- 
tionary models  to  a sequence  of  such  intervals.  The  correspond- 
ing spectra  are  completely  described  by  the  AR  coefficients 
and  the  variance  a 2 , and  are  given  by 

2o  2 

P(f)  = 0<f  < ^ (2) 

D — — 

I 1-  1 a,  exp(-2TTifk)  j 2 
k=l  K 

The  order  p=p(t)  of  the  model  has  been  estimated  using  the 

Final  Prediction  Error  criterion  and  is  a function  of  time. 

The  time  variation  of  the  spectrum  is  now  determined  by  the 

behavior  in  time  of  a.(t),...,a  . (t)  and  o2 (t) . Usuallv 

i p(t)  z 

the  time  variation  is  described  in  terms  of  selected  snectral 
peaks,  but  this  leads  to  loss  of  information. 

The  autoregressive  modelling  was  used  to  examine  the  changes 
in  the  LP-spectrum  during  two  microseismic  storms  described 
by  Korhonen  and  Pirhonen  (1976).  The  first  example  is  a cvclone 
center  over  the  North  Sea  in  the  time  period  Parch  1°72. 

An  AR-model  of  order  7 produced  a satisfactory  approximation 
to  the  spectrum.  However,  most  of  the  spectral  information  is 
contained  in  the  first  two  coefficients.  The  white  noise 
variance  acts  as  a scaling  factor.  The  time  development 


«' 
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of  the  coefficient  a,  and  the  variance  o for  the  vertical 

1 z 

component  of  the  LP-noise  is  shown  in  Fig.  VII. 8.1  a)  and  h) 
respectively.  The  parameter  has  a sharo  peak  on  March  23 
at  about  18h00  for  all  three  components.  This  corresponds 
to  the  maximum  power  of  the  first  and  second  peak  (Korhonen 
and  Pirhonen,  1976)  of  the  vertical  component  spectrum. 

It  is  also  interesting  to  note  that  the  wind  velocity  has 
a maximum  between  March  22  and  March  23,  which  corresponds 
to  the  time  interval  of  greatest  change  for  the  coefficient 
a..  For  more  details  and  for  a discussion  of  the  second 
example  (a  cyclone  center  approaching  the  North  Atlantic 
in  the  time  interval  May  24-27  1972)  we  refer  to  Sandvin 
and  Tjdstheim  (1977). 


O.A.  Sandvin 
D.  Tjpstheim 
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VII. 9 A Pattern  Recognition  Approach  to  Seismic  Discrimination. 


Part  II:  Classif ication 

As  explained  in  Part  I (Tj0stheim,  1976)  the  seismic  discrimina- 
tion problem,  when  identified  as  a pattern  recognition  problem, 
can  be  separated  into  two  stages:  feature  extraction  and  clas- 
sification. The  construction  of  a primary  feature  vector  Y 
containing  short  and  long  period  features  was  described  in 
Tj0stheim  (1976).  Using  principal  component  analysis  the  vector 
Y has  been  further  reduced  to  an  8-dimensional  feature  vector 
Z where 


2 = lZEX(li 


. , ZEX ( 4 ) , ZEQ( 1) , 


' ' ZEQ  ^ lj 


(1) 


with 


1 ’ 


ZEX(l)  - ' — i , EX 


(2) 


where  ft.  „ , i=l,...,4  are  the  four  major  estimated  principal 

l f IliA 

component  vectors  of  the  explosion  data  set.  The  quantities 
ZEQ(i)f  i=l,...,4  are  defined  similarly. 


It  is  assumed  that  the  Z vector  is  governed  by  a multivariate 
Gaussian  probability  density  function  P^(z)  or  P2(z)  cor- 
responding to  the  explosion  (1)  population  or  earthquake  (2) 
population.  The  density  function  estimation  problem  then 
reduces  to  estimating  the  respective  mean  vectors  and  co- 
variance  matrices.  The  following  decision  rule  was  applied: 

For  a given  observed  event  with  a corresponding  feature  vector 
Z,  the  quantities  P^(j2)  and  P2  (Z)  were  computed,  and  Z was 
assigned  to  the  population  i having  the  largest  probability 
P • ( 2, ) of  having  Z_  occurring. 

In  practice  this  means  that  the  decision  rule  is  to  declare 
an  explosion  if  S^Z)  - §2(Z)  is  positive  and  an  earthquake 
if  this  difference  is  negative.  Here 
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s.(z)  = log!fti;  -^(z-u±)  1(z-0±)T  (3) 

where  and  (m,  i=l,2  are  the  explosion  and  earthquake 
Z-vector  covariance  matrices  and  mean  vectors  respectively. 

This  classification  rule  was  applied  on  the  data  set  of  52 
explosions  and  73  earthquakes  described  in  Tj0stheim  (1976).  The 
results  are  displayed  in  the  bottom  part  of  Fig.  VII .9.1  where 
the  events  have  been  plotted  according  to  their  value  of 
^l(*)~^2^*^*  i-s  seen  that  all  events  are  correctly  classified. 

For  matters  of  comparison  we  treated  the  m.  :M  data  and  the 

D S 

Xl:X2  discriminant  data  of  Tj0stheim  and  Husebye  (1976)  in 
a similar  way.  That  is,  we  used  the  discriminant  (3)  but  with 
Z replaced  by  U = (m^M  ) or  V = (XI,  X2)  . (The  notation  Ll:L2 
is  used  on  the  figure.)  From  Fig.  VII. 9.1  it  is  seen  that  the 
mb:Ms  ^ata  Produces  serious  overlap  and  that  the  performance 
is  clearly  inferior  to  the  discriminant  based  on  § ( Z)  0 ( Z)  . 

The  Xl:X2  data  produces  complete  separation,  but  with  one 
explosion  very  close  to  the  explosion-earthquake  decision  line. 

D.  Tj0stheim 
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VII. 10  A New  Method  of  Spectral  Estimation  for  Spatial  Data 

Geophysical  data  in  general  and  seismolog ical  data  in  particular 
are  often  recorded  both  in  space  and  in  time.  While  a number 
of  efficient  statistical  techniques  exist  for  analysis  of 
data  observed  in  time,  this  is  not  so  for  spatial  data.  We 
propose  here  a new  type  of  analysis  for  spatial  variables 
observed  on  a grid.  The  proposed  technique  might  be  viewed 
as  a generalization  of  the  time-series  autoregressive  analysis, 
but  there  is  also  some  justification  for  considering  it  a 
spatial  extension  of  Burg's  (1975)  maximum  entropy  analysis. 

For  simplicity  consider  the  case  where  we  have  a geophysical 
quantity  F(x^,x2)  observed  on  a regular  grid  in  the  X^-X2 
plane.  We  use  the  notation  S<x,y]  for  the  vectors  of  integers 
u=(u^,u2)  which  are  such  that  xk<u^<yk  for  k=l,2  but  u*x,  and 
seek  to  model  F(x)  = F(x^,x2)  by  a one-sided  autoregressive 
model  in  the  plane 

F (x)  - l a (y)  F(x-y)  = Z (x)  (1) 

y€S<0 , p] 

where  Z(x)  = Z(x^,x2)  is  a spatial  white  noise  series.  Cor- 
responding to  the  time  series  case  we  obtain  an  estimate 
£(A)  = £(>^,X2)  of  the  two-dimensional  power  spectral  density 
of  F ( x ) as 


f ( > ) = 


jl  - /.  a (y)  exp{  -i  [ y , x ] } | 2 

y£S<0 , p] 


We  have  done  some  preliminary  simulation  experiments,  using 
artificially  generated  data  F(x).  In  each  case  we  computed 
f(A)  as  well  as  an  FFT  power  spectral  estimate.  Fig.  VII.  10.1 
shows  the  results  for  a model  of  the  form 


F(Xi,x2)  Z(x1,x2)+0.85o2  cos(1.5x1+1.5x2)+0.85oz  cos  1.5x2  (3) 
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where  we  have  tv/o  cosines  embedded  in  two-dimensional  spatial 
white  noise.  It  is  seen  that  the  two  spectral  peaks  are  much 
sharper  in  the  autoregressive  case.  This  suggests  that  as  in 
the  time  series  case  there  may  be  cases  where  spatial  auto- 
regressive procedures  are  superior  to  FFT  estimates.  More 
tests  are  required  to  put  this  conclusion  on  a firmer  basis. 

D.  Tjpstheim 
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Estimating  the  Autocorrelation  Function  of  Spatio-Temporal 
Variables 

In  geological  exploration  as  well  as  in  air  pollution  studies 
it  is  often  necessary  to  interpolate  from  a sparse  and  irregu- 
larly positioned  set  of  measurement  points.  The  principal 
practical  difficulty  (Delfiner,  1976)  in  studies  of  this  sort 
consists  in  obtaining  good  estimates  for  the  spatial  auto- 
correlation function.  Motivated  by  these  problems  we  suggest 
(Tj0stheim,  1977)  an  estimate  for  the  autocorrelation  function 
fo<r  spatio-temporal  variables.  The  estimate  is  meaningful  for 
non-integer  lags  and  for  irregularly  positioned  data  points. 

It  is  constructed  using  a double  Fourier  transform  technique. 
The  statistical  properties  of  the  estimate  are  studied  in  a 
special  case.  For  more  details  and  for  a brief  summary  of 
applications  to  air  pollution  modelling,  we  refer  to  Tj0stheim 
(1977)  . 

D.  Tj0stheim 
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